BDF-α: A Multistep Method with Numerical Damping Control
نویسندگان
چکیده
When solving numerically the stiff second order ODE system obtained after semidiscretizing the wave-type partial differential equation (PDE) with the finite element method (FEM), and similarly to the HHTα method, which allows the numerical damping of the undesirable high frequency modes associated to FEM semidiscretization, we have constructed a modification of the 2-order BDF method (the BDF2 method), which we have called BDF-α. This new method is second-order accurate and with a smaller local truncation error than the BDF2, it is unconditionally stable for some values of α and it permits a parametric control of numerical dissipation.
منابع مشابه
The integration of stiff systems of ODEs using multistep methods
The second order Ordinary Differential Equation (ODE) system obtained after semidiscretizing the wavetype Partial Differential Equation (PDE) with the Finite Element Method (FEM), shows strong numerical stiffness. Although it can be integrated using Matlab ode-solvers, the function ode15s offered by Matlab for solving stiff ODE systems does not result very efficient as its resolution requires t...
متن کاملA-Stable Time Discretizations Preserve Maximal Parabolic Regularity
It is shown that for a parabolic problem with maximal Lp-regularity (for 1 < p < ∞), the time discretization by a linear multistep method or Runge–Kutta method has maximal `p-regularity uniformly in the stepsize if the method is A-stable (and satisfies minor additional conditions). In particular, the implicit Euler method, the Crank–Nicolson method, the second-order backward difference formula ...
متن کامل1 Week 11 : 4 and 6 November 2008
We now discuss the theory of fixed–order constant–step–size linear multistep methods. This classical theory is presented in [HNW93] Sect. III.1–4 (nons-tiff) and [HW96] Sect. V.1 (stiff), as well as [Gea71] Chap. 7, [Lam91] Chap. 3, [Sha94] Chap. 4, and many other sources. While multistep methods are rarely used in practice with fixed order and constant step size, this theory is useful because ...
متن کاملA-EBDF: an adaptive method for numerical solution of stiff systems of ODEs
In this paper a one parameter predictor–corrector method, which we call it A-EBDF, is introduced and analyzed. With a modification of A-BDF and EBDF methods we propose a multistep method whose region of absolute stability is larger than those of A-BDF and EBDF methods. © 2004 IMACS. Published by Elsevier B.V. All rights reserved. MSC: 65L05
متن کاملMultistep Methods Using Higher Derivatives and Damping at Infinity*
Linear multistep methods using higher derivatives are discussed. The order of damping at infinity which measures the stability behavior of a fc-step method for large h is introduced. ,4-stable methods with positive damping order are most suitable for stiff problems. A method for computing the damping order is given. Necessary and sufficient conditions for A -stability, A (a) -stability and stif...
متن کامل